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ABSTRACT 

Context. How complex organic - and potentially prebiotic - molecules are formed in regions of low- and high-mass star-formation 
remains a central question in astrochemistry. In particular, with just a few sources studied in detail, it is unclear what role environment 
plays in complex molecule formation. In this light, a comparison of relative abundances of related species between sources might be 
useful to explain observed differences. 

Aims. We seek to measure the relative abundance between three important complex organic molecules, ethylene glycol((CH 20 H) 2 ), 
glycolaldehyde (CH 2 OHCHO) and methyl formate (HCOOCH 3 ), toward high-mass protostars and thereby provide additional con¬ 
straints on their formation pathways. 

Methods. We use IRAM 30-m single dish observations of the three species toward two high-mass star-forming regions - W51/e2 and 
G34.3-I-0.2 - and report a tentative detection of (CH 20 H )2 toward both sources. 

Results. Assuming that (CH 20 H) 2 , CH 2 OHCHO and HCOOCH 3 spatially coexist, relative abundance ratios, HC 00 CH 3 /(CH 20 H) 2 , 
of 31 and 35 are derived for G34.3-I-0.2 and W5 l/e2, respectively. CH 2 OHCHO is not detected, but the data provide lower limits to the 
HCOOCH 3 /CH 2 OHCHO abundance ratios of >193 for G34.3-I-0.2 and >550 for W51/e2. A comparison of these results to measure¬ 
ments from various sources in the literature indicates that the source luminosities may be correlated with the HC 00 CH 3 /(CH 20 H )2 
and HCOOCH 3 /CH 2 OHCHO ratios. This apparent correlation may be a consequence of the relative timescales each source spend 
at different temperatures-ranges in their evolution. Furthermore, we obtain lower limits to the ratio of (CH 20 H) 2 /CH 20 HCH 0 for 
G34.3-I-0.2 (> 6 ) and W51/e2 (> 16). This result confirms that a high (CH 20 H) 2 /CH 20 HCH 0 abundance ratio is not a specific 
property of comets, as previously speculated. 

Key words, astrochemistry - ISM: molecules - ISM: abundances - ISM: individual objects: W51/e2, G34.3-I-0.2 - 
line: identification - methods: observational astrochemistry 


1. Introduction 

A central question in astrochemistry is how, where, and when 
complex organic molecules form. Although the total number of 
detected molecules in the interstellar medium^ (ISM) continu¬ 
ously increase (e.g. Herbst & van Dishoeck 2009), the answers 
to the question above still remain unclear. To date there is no con¬ 
sensus on how complex organic molecules form in dense regions 
of the ISM, despite the increasing number of detections. One 
promising suggestion was that warm gas-phase chemistry, fol¬ 
lowing evaporation of simple ices, could be a primary formation 
pathway (e.g. Millar et al. 1991; Charnley et al. 1992). However, 
more recent laboratory experiments and chemical modeling has 
shown that this mechanism is not effective enough to account for 
the observed abundances (e.g. Geppert et al. 2006). An alterna¬ 
tive formation mechanism involves UV induces radicals. Garrod 
et al. (2008) propose that radicals, during the warm up phase, can 
migrate on the grains surfaces and form complex species which 
are then released into the gas-phase at higher temperatures. The 
initial ice composition and the amount of UV radiation have 
also proved to play an important part of the formation process 
(Oberg et al. 2009). In order to determine the formation path¬ 
way of COMs, it is useful to determine abundance ratios between 


http://www.astro.uni-koeln.de/cdms/molecules 


different related species as these ratios, in comparison with the 
predicted abundance ratios from chemical models, can provide 
constraints on the formation processes. Indeed, variations in the 
abundance profiles can reflect the physical and chemical condi¬ 
tions that are occurring. It is therefore important to observe these 
species in different environments. 

Some of the simplest species in this context include the 
oxygen-bearing complex organic molecules associated to gly¬ 
colaldehyde, CH 2 OHCHO, including its isomer methyl formate 
(HCOOCH 3 or CH 3 OCHO) and the reduced alcohol version 
of CH 2 OHCHO, ethylene glycol, (CH 20 H )2 (also commonly 
known as anti-freeze). By constraining the relative abundances 
of these species in different environments the hope is to be able 
to explore, e.g., the importance of initial chemical conditions, 
temperature and irradiation in the formation for comparison, 
e.g., to laboratory experiments (e.g. Oberg et al. 2009) and as in¬ 
put for sophisticated chemical models (e.g., Garrod et al. 2008). 
For example, in their laboratory experiments Oberg et al. (2009) 
show that the relative abundances of HCOOCH 3 to CH 2 OHCHO 
and (CH 20 H )2 are strongly dependent on both the ice tempera¬ 
ture and exact ice composition in terms of the relative amounts 
of CO andCH 30 H. 

So far, (CH 20 H )2 has been detected toward high-mass 
sources such as the Galactic center source Sgr B2(N) by Hollis 
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Table 1: Transitions of (CH 20 H )2 in the observed frequency range. 


Transition 

Frequency 

E„p 


Rotational diagram 


[MHz] 

[K] 

[s-‘] 


3 mm observations 

o 

II 

oo 

1 

II 

6 \ 

100333.6424 

23.6 

-4.6773 

W51/e2 

93 6 V — 1 — 83 5 V — 0 

100490.6121 

27.0 

-4.6800 


111.10 V = 1 — 102,9 V = 1 

102064.0770 

33.9 

-4.6692 


lOl.io V = 1 - 9 i,9 V = 0 

102291.5860 

26.2 

-4.5798 


92,7 V = 1 — 82,6 V = 0 

102539.4267 

25.0 

-4.6355 


lOo.io V = 1 - 9 o,9 V = 0 

102689.8384 

26.1 

-4.4419 


119,2 V = 0 - 109,1 V = 1 

105574.9490 

72.2 

-4.9342 


118,4/3 V = 0 — 108,3/2 V = 1 

105644.1940 

63.9 

-4.5179 


111,10 V = 0 - 10i,9 V = 1 

105701.0020 

33.6 

-4.3644 


117,5/4 V = 0 — 107,4/3 V = 1 

105735.3520 

56.5 

-4.4051 


113,9 V = 0 — 103,8 V = 1 

105834.6960 

36.9 

-4.5676 


116,6/5 V = 0 - 106,5/4 V = 1 

105866.9546 

50.1 

-4.3228 


115,7 V = 0 — 105,6 V = 1 

106072.0524 

44.7 

-4.6197 


115,6 V = 0 - 105,5 V = 1 

106088.9178 

44.8 

-4.5104 


1 mm observations 

24i7,8/7 V = 0 - 2317,7/6 V = 1 

238828.3252 

289.3 

-3.4806 


24i6,9/8 V = 0 - 2316,8/7 V = 1 

238904.3833 

273.1 

-3.4330 


24i5,io/9 V = 0 - 23i5,9/8 V = 1 

238993.6343 

257.9 

-3.3925 


30i2,19/18 V = 1 - 30ii,19/20 V = 0 

238994.0910 

299.7 

-4.8643 


225,17 V = 1 — 215,16 V = 0 

238994.7545 

138.6 

-3.5149 


24i3,i2/11 V = 0 - 23i3,ii/io V = 1 

239231.3329 

230.4 

-3.3269 


24i2,13/12 V = 0 - 23i2,12/11 V = 1 

239395.0677 

218.2 

-3.3002 


24ii,14/13 V = 0 - 23ii,13/12 V = 1 

239605.3330 

206.9 

-3.2765 


25 3,23 V = 0 — 243,22 V = 1 

239792.7984 

161.6 

-3.4697 


24io,15/14 V = 0 - 23io,14/I3 V = 1 

239883.5631 

196.7 

-3.2555 

W51/e2 

252,23 V = 0 — 242,22 V = 1 

239957.1582 

161.6 

-3.4917 


234,20 V = 1 — 224,19 V = 0 

239980.0390 

144.1 

-3.4495 


243,21 V = 0 — 233,20 V = 1 

240147.9627 

155.2 

-3.4704 


249,16/15 V = 0 - 239,15/14 V = 1 

240265.7687 

187.4 

-3.2365 


224,19 V = 0 — 213,18 V = 0 

240608.1463 

132.6 

-3.9839 


251/0,25 V = 1 - 24i/o,24 V = 0 

240778.3210 

148.0 

-3.1643 


248,17 V = 0 - 238,16 V = 1 

240807.8795 

179.2 

-3.5201 

W51/e2and G34.3+0.2 

248,16 V = 0 - 238,15 V = 1 

240828.8863 

179.2 

-3.5199 


224,18 V = 1 — 214,17 V = 0 

240875.3680 

135.5 

-3.4905 


59io,5o V = 1 - 599,51 V = 1 

240888.3533 

925.5 

-4.4992 


245,20 V = 0235,19 V = 1 

241291.2695 

160.7 

-3.5046 


247,18 V = 0 - 237,17 V = 1 

241545.2626 

172.1 

-3.5045 


247,17 V = 0 — 237,16 V = 1 

241817.6908 

172.1 

-3.5026 


246,19 V = 0 - 236,18 V = 1 

241860.7330 

166.0 

-3.5459 


23i5,9/8 V = 1 - 22i5,8/7 V = 0 

242244.6928 

246.4 

-3.4004 


23i4,io/9 V = 1 — 22i4,9/8 V = 0 

242246.3392 

232.2 

-3.3608 


23i6,8/7 V = 1 - 22i6,7/6 V = 0 

242267.1576 

261.6 

-3.4470 


23i3,ii/io V = 1 - 22i3,io/9 V = 0 

242277.7167 

218.9 

-3.3268 


23i7,7/6 V = 1 - 22i7,6/5 V = 0 

242309.7823 

277.8 

-3.5027 


23i2,12/11 V = 1 - 22i2,11/10 V = 0 

242347.0348 

206.7 

-3.2973 


23i8,6/5 V = 1 - 22i8,5/4 V = 0 

242369.7298 

295.0 

-3.5709 


262/1,25 V = 0 - 252/1,24 V = 1 

242389.2838 

166.9 

-3.1601 


23ii, 13/12 V = 1 - 22ii,12/11 V = 0 

242466.6909 

195.4 

-3.2714 


23io,14/13 V = 1 - 22io,13/12 V = 0 

242656.2283 

185.2 

-3.2486 

W51/e2and G34.3+0.2 

233,20 l-22319v = 0 

242897.4293 

143.6 

-3.7464 


239,15/14 V = 1 - 229,14/13 V = 0 

242948.2909 

175.9 

-3.2282 

W51/e2 

246,18 V = 0 - 236,17 V = 1 

243259.7398 

166.3 

-3.8001 


238,16 V = 1 — 228,15 V = 0 

243396.0258 

167.7 

-3.5106 


238,15 V = 1 — 228,14 V = 0 

243408.4357 

167.7 

-3.5105 


235,19 V = 1 — 225j8 V = 0 

243636.5668 

149.1 

-3.4704 


237,17 V = 1 — 227 16V = 0 

244054.2878 

160.5 

-3.4934 


237,16 V = 1 — 227,15 V = 0 

244233.4545 

160.5 

-3.4927 


243,22 V = 1 — 233,21 V = 0 

244399.6722 

150.1 

-3.4754 


242,22 V = 1 — 232,21 V = 0 

244685.1477 

150.1 

-3.4399 


236,18 V = 1 — 226,17 V = 0 

244879.9193 

154.4 

-3.5188 


27i/o,27 V = 0 - 26i/0,26 V = 1 

245022.7624 

171.4 

-3.1407 


246,19 V = 1 - 236,17 V = 1 

245410.9218 

166.3 

-3.7580 



Notes. Catalog values for the (CH 20 H )2 transitions that can be excited in the observed frequency range for W51/e2 and G34.3+0.2. Some faint 
transitions have been excluded from this table, please section 3. The source names are listed for the lines that are used in the rotational diagram 
analysis. 
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et al. (2002), marginally toward W51 el/e2 by Kalenskii & Jo¬ 
hansson (2010) and recently also toward the Orion Kleinmann- 
Low nebula by Brouillet et al. (2015). (CH 20 H )2 has also been 
observed toward the low-mass Class 0 protostars IRAS 16293- 
2422 (Jprgensen et al. 2012, in prep.) and NGC 1333 IRAS 
2A (Maury et al. 2014; Coutens et al. 2015) as well as toward 
the intermediate-mass protostar NGC 7129 FIRS2 (Fuente et al. 
2014). CH 2 OHCHO has been detected toward Sgr B2(N) (Hollis 
et al. 2000, 2001, 2004; Halfen et al. 2006; Requena-Torres et al. 
2008), the high-mass hot molecular core G31.41-1-0.31 (Beltran 
et al. 2009), in IRAS 16293-2422 by Jprgensen et al. (2012) 
and recently also in NGC 7129 FIRS2 (Fuente et al. 2014), 
NGC 1333 IRAS 2A (Coutens et al. 2015; Taquet et al. 2015) 
and NGC 1333 IRAS 4A (Taquet et al. 2015). HCOOCH 3 is 
the most abundant isomer of CH 2 OHCHO and has previously 
been observed in numerous hot cores and corinos, (e.g. Blake 
et al. 1987; Bisschop et al. 2007; Demyk et al. 2008; Favre 
et al. 2014, 2011). There are some notable differences in terms 
of the abundance ratios. For example, Coutens et al. (2015) find 
a (CH 20 H) 2 /CH 20 HCH 0 ratio of ~5 in NGC 1333 IRAS 2A, 
while Jprgensen et al. (2012 and in prep.) find a lower value of 
~ 1 in IRAS 16293-2422. These changes hint that it might be use¬ 
ful to explore these ratios in different sources sampling similar 
physical conditions. 

This paper presents IRAM 30-m observations of (CH 20 H )2 
and CH 2 OHCHO toward the high mass protostars W51/e2 (dis¬ 
tance = 5.4 kpc, Sato et al. 2010) and G34.3H-0.2 (distance = 
3.8 kpc, Kurtz et al. 2000). The aim of this study is to deter¬ 
mine the relative abundance of (CH 20 H )2 and CH 2 OHCHO to 
HCOOCH 3 and compare to other sources from the literature. 
Overall our aim is to investigate the use of abundance ratios of 
species that are believed to be chemically related to explore the 
origin of complex molecules in the dense interstellar medium. 
In Section 2, the IRAM 30-m observations are presented. Data 
modeling, results and analysis are presented in Section 3 and 
discussed in Section 4. 


2. Observations 

The observations were performed with the IRAM 30-m tele¬ 
scope at Pico Veleta, Spain on December 13 and 14, 2012 
for W51/e2 and G34.3-I-0.2, respectively. The coordinates of 
the phase tracking center used for the two sources were ay 2000 
= 19''23'”43?9, 5/2000 = 14°30'34'.'8 for W51/e2 and ay 2000 
= 18'‘53'"18!6, 5/2000 = 0ri4'58'.'0 for G34.3H-0.2. The ob¬ 
servations were performed in position switching mode using 
[-600", 0"] as reference for the OFF positions. The spectral 
setup was chosen to target (CH 20 H )2 and CH 2 OHCHO in the 
observations. The EMIR receiver was used in dual band polar¬ 
isation (E090/E230) in connection with i) the 200 kHz Eourier 
transform spectrometer (ETS) back-end in the frequency ranges 
99.70-106.3 GHz and 238.2-246.0 GHz for the E090 and E230 
bands, respectively; ii) the WILMA back-end in the frequency 
ranges 99.88-103.6 GHz and 238.4-242.1 GHz for the E090 
and E230 bands, respectively. Data reduction was performed us¬ 
ing The Continuum and Line Analysis Single-dish Software^ 
(CLASS). The resulting ETS spectra were smoothed to a spec¬ 
tral resolution of 1.15 km s ' for the E090 band and 1.22 km s * 
for the E230 band while the WILMA spectra were smoothed to 
5.88 km s ' and 2.49 km s ' for the E090 and E230 bands, re¬ 
spectively. Spectra resulting from the ETS back-end presented a 
standing wave pattern and the WILMA observations were there- 

2 http://www.iram.fr/IRAMFR/GILDAS 


fore used as a sanity check for the LTS data and to confirm the 
detections made in the LTS observations for the lines where the 
frequency of the two sets of observations match. The standing 
wave present in the LTS observations was removed using a fast 
Fourier transform in the data reduction. 

Throughout this paper, the intensity is given as the main 
beam brightness temperature (r^b), which is defined as 


where Ta is the antenna temperature, Egfr is the forward 
efficiency and Beft is the beam efficiency. The values used 
were Leff = 0.94 and 0.92 and Beft = 0.78 and 0.58 for 
1 mm and 3 mm respectively. The half-power beam sizes are 
~ 10" and ~ 24" for the observations at 1 mm and 3 mm, 
respectively. The reduced data-sets are available for download at 
http://youngstars.nbi.dk/projects/HighMass_organics/index.html. 

3. Analysis and Results 

Spectra from both sources show a rich forest of lines character¬ 
istic for high-mass sources. A total of 21 and 19 lines have peak 
temperature above 5 K for W51/e2 and G34.3H-0.2, respectively; 
the strongest lines in both spectra is the CS 5^ transition at 
244935 MHz. It is important to know that there is still a remnant 
oscillation in the baseline even after the removal of the standing 
wave in the spectra. This increases the RMS noise of the spec¬ 
tra, which complicates the analysis of faint lines. Therefore, in 
addition to the global baseline subtraction, an additional local 
zero or first order baseline subtraction was therefore performed 
before making a Gaussian fit of the (CH 20 H )2 lines. This addi¬ 
tional baseline subtraction is applied in a range of + 100 km s ' 
for each line. As for the data reduction, CLASS was used for the 
additional baseline subtraction as well as Gaussian profile fits 
to the detected lines. The resulting RMS of the local baseline is 
~ 60 mK for the 1 mm observations and ~ 7 mK for the 3 mm 
observations. 

In order to ensure proper line identification, we have checked 
the line observations against entries in Splatalogue database for 
astronomical spectroscopy^. In addition, we have made a ref¬ 
erence model where we produced synthetic spectra for the line 
emission of common species (CH 3 OCH 3 , HCOOH, CH 3 CHO, 
CH 3 OH, C 2 H 5 OH and CH 3 CN) in order to visually exclude 
lines that are blended with any of these molecules. Table 1 lists 
spectroscopic parameters for the (CH 20 H )2 transitions that can 
be excited in the observed frequency range. Transitions with 
log(Aui) < -5 for the 3 mm observations and log(Aui) < -4 and 
Bup > 300 K for the 1 mm observations have been excluded from 
the table, as these transitions are predicted from the synthetic 
spectra to have peak intensities ^ 0.02 K, i.e., not detectable. 

The spectroscopic data for (CH 20 H )2 come from Christen 
& Muller (2003) and Christen et al. (1995) and are available 
from the CDMS database"' (Muller et al. 2001, 2005), while the 
spectroscopic data for HCOOCH 3 and CH 2 OHCHO are from 
Ilyushin et al. (2009) and Carroll et al. (2010) respectively, avail¬ 
able from the JPL database^ (Pickett et al. 1998). For the anal¬ 
ysis, the (CH 20 H )2 lines that are reasonably well separated and 
have a peak temperature above 3 cr were selected. For W51/e2 
we have included the lines at ~ 100333 and ~ 242656 MHz al¬ 
though they are partially blended, as it is possible to distinguish 

3 www.splatalogue.net 
* http://www.astro.uni-koeln.de/cdms 
^ http://spec.jpl.nasa.gov/home.html 
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Table 2: Transitions of (CH 20 H )2 observed toward W5 l/e2 and G34.3+0.2 with the IRAM 30-m telescope. 


Transition 

Freq. 

E„p 

logioA„/ 

J TmbAv 

^ LSR 

Av 

Peak T„fc 

4,.k,. V" 

[MHz] 

[K] 

[s-'l 

[ K km s“'] 

[km s“'] 

[km s“'] 

[K] 

W51/e2 

00 

II 

1 

00 

kj 

< 

II 

0 

100333.64 

24 

-4.6773 

0.2 

55.1 

5.5 

0.03 

24io, 15/14 V = 0 - 23l0, 14/13 V = 1 

239883.56 

197 

-3.2555 

2.7 

55.6 

6.8 

0.38 

248,17 V = 0 - 238,16 V = 1 

240807.88 

179 

-3.5201 

2.0 

56.7 

5.9 

0.32 

23io,14/13 V = 1 - 22io,13/12 V = 0 

242656.23 

185 

-3.5496 

4.0 

56.9 

7.4 

0.50 

239 , 15/14 V = 1 - 229 , 14/13 V = 0 

242948.29 

176 

-3.5292 

1.4 

56.5 

5.0 

0.27 

G34.3+0.2 

248,17 V = 0 - 238,16 V = 1 

240807.88 

179 

-3.5201 

1.0 

57.7 

4.1 

0.23 

23io,14/13 V = 1 - 22io,13/12 V = 0 

242656.23 

185 

-3.5496 

1.0 

58.4 

5.0 

0.19 


Notes. The table lists spectroscopic parameters and observed quantities from Gaussian fits to the detected (CH 20 H )2 lines: integrated line intensity 
(f TmbAv), line position (Ylsr), line width (Av) and peak temperature (Peak T^i,). The error of J T^bAv is estimated to ~ 30% as determined by 
the observational uncertainty, while the errors on Ylsr and Av are ~ 1.0 and ~ 2.0 km s“', respectively. The line intensities have been added 
together for the hyper fine structure transitions (denoted by / in the quantum numbers). 


the (CH 20 H )2 peak from the other peaks and perform a multiple 
Gaussian fit which includes all relevant peaks. 

Figure 1 and 2 show a zoom-in of the area of ~ 200 km s“\ 
which corresponds to ~ 160 MHz in the 1 mm observations 
and ~ 70 MHz in the 3 mm observations, around each detected 
line in both sources after the local baseline subtraction. Super¬ 
imposed onto the observed spectra are the synthetic spectra of 
(CH 20 H )2 as well as the other species investigated here, as to 
demonstrate that the chosen (CH 20 H )2 lines are not blended. 
The fit to the entire observed spectra for both sources are shown 
in Figure B.1-B.8 in the appendix. Table A.l in the appendix lists 
the estimated column densities of the molecules in the reference 
model. Table 2 lists spectroscopic parameters and the observed 
quantities from the fits of the (CH 20 H )2 lines: the integrated line 
intensities (J TnibAv), position (Vlsr), line width (Av) and peak 
temperature (Peak Tnib)- 

A total of 35 and 52 well separated lines have been detected 
above the 3 cr level in the 1 mm and 3 mm data for HCOOCH 3 
in W51/e2 and G34.3+0.2 respectively. For (CH 20 H) 2 , 4 and 2 
lines, toward W51/e2 and G34.3+0.2 respectively, are detected 
in the 1 mm observations; while only one line is detected at 
3 mm toward for W5 l/e2. Many other (CH 20 H )2 lines as well as 
some potential CH 2 OHCHO lines are present in both data-sets, 
but they are either blended with other species or too faint to be 
properly detected above a 3 cr limit. 

3.1. Modeling 

Assuming LTE and optically thin emission, the rotational dia¬ 
gram method can be used to determine the column density, Atot, 
and the rotational temperature, T^ot (Goldsmith & Ganger 1999). 
The approach used in this work follows the formalism described 
by Goldsmith & Ganger (1999) with the following assumptions: 

1. (CH 20 H) 2 , CH 2 OHCHO and HCOOCH 3 are emitting from 
the same region and thus have the same source size, ^source- 

2. (CH 20 H) 2 , CH 2 OHCHO and HCOOCH 3 are in LTG, which 
implies that the excitation temperature, Tex, is equal to the 
kinematic temperature ^kin. for all three species. Using the 
rotational diagram method one obtains the rotational tem¬ 
perature which in LTG is Tkin = T^ot - Tex- 

3. The observations at 1 mm and 3 mm trace the same gas. 

According to our assumptions, the source size and Trot derived 
from the HCOOCH 3 analysis can be used to derive the column 


density of (CH 20 H )2 and set an upper limit on the column den¬ 
sity of CH 2 OHCHO. A possible caveat is that the large beam size 
of our observations compared to the source sizes is not guaran¬ 
teeing that the measured lines arise from species present in the 
same gas. Specifically, we obtain similar mean values of Vlsr 
for (CH 20 H )2 and HCOOCH 3 in both sources (56.4 + 0.5 km s ' 
and 56.0+0.8 km s“' in W51/e2 along with 57.7+0.5 km s“' and 

58.4 + 1.2 km s ' in G34.3+0.2 for (CH 20 H )2 and HCOOCH 3 , 
respectively) so here the assumption appears reasonable. 

The same formalism in Goldsmith & Ganger (1999) used 
to create a rotational diagram can also the used to generate a 
synthetic spectrum of the line emission of a specific molecules 
by using the source size, ^source, the rotational temperature. Trot, 
the line width, Av, and the column density, Atot, as input pa¬ 
rameters. It is possible to correct for the optical depth of the 
line emission in case it deviates from being completely optically 
thin (Goldsmith & Ganger 1999). A synthetic spectrum gener¬ 
ated with inputs parameters derived from the rotational diagram 
therefore serves as a check for self-consistency of the result. We 
use adopted source sizes from the literature and determine the ro¬ 
tational temperature and column densities in the following anal¬ 
ysis. 

The observed width of the detected lines have been deter¬ 
mined from the Gaussian fit. Gor HCOOCH 3 in W51/e2 we ob¬ 
tain mean values of Av = 6.4 + 0.8 km s“' in the 1 mm data and 
Av = 6.1 + 1.1 km s“' in the 3 mm data. In the same source, the 
mean value of the line width of the (CH 20 H )2 lines is 6.3 + 0.9 
km s ' for the 1 mm data and the only detected value is 5.5 km 
s“' in the 3 mm data. Gor G34.3+0.2 we find Av = 6.5 + 1.5 km 
s“' in the 1 mm data and Av = 7.3 + 1.3 km s“' in the 3 mm data 
for HCOOCH 3 and Av = 4.5 + 0.5 km s ' for the (CH 20 H )2 
lines in the 1 mm data. In the below analysis, a fixed line width 
of 6.0 km s ' was chosen as input parameter in synthetic spectra 
for all species. 

3.2. HCOOCH 3 

In the rotational diagram analysis for W51/e2, a source size of 

2.4 ” (Zhang et al. 1998) is applied and results in a rotation tem¬ 
perature of Trot = 120 K and a column density of 1.1 x 10' * cm“^. 
We estimate a 40 K uncertainty on the rotational temperature, 
which contribute to a ~ 20% error on the column density. This, 
combined with an observational uncertainty of ~ 30%, returns 
an estimated overall error of 30 ^0 % for the column density. 
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Fig. 1 : Observed (black) and synthetic (red) spectra of the five transitions of (CH20H)2 in W51/e2 which are used for the analysis. Also shown 
are the synthetic spectra of all other species investigated in this study (blue) to demonstrate that the (CH 20 H )2 lines in question are not blended. 
The transitions are located at 239883.56, 240807.88, 242656.23, 242948.29 and 100333.64 MHz - marked by black vertical lines in the plot. The 
displayed frequency region in each plot corresponds to ~ 200 km s^F 


Table 3: Input parameters for synthetic spectra 


Parameter 

G34.3+0.2 

W51/e2 

^source [ ] 

7.6 

2.4 

Tro,[K] 

140 

120 

Av [kms^'] 

6.0 

6.0 

HCOOCH3 Ntot[cm- 2 ] 

5.8 X 10'® 

1.1 X 10'^ 

(CH 20 H )2 (V,ot[cm-2] 

1.9 X 10'® 

3.1 X 10 '® 

CH2OHCHO Ntot[cm 2 ] 

< 3 X 10'4 

< 2 x 10 '® 


Notes. The source sizes, rotational temperatures, line widths and col¬ 
umn densities which are used as input parameters when creating syn¬ 
thetic spectra. As described in the text, the estimated uncertainties are 
±40 K for Trot and ±30-40% for the column density. 


For G34.3+0.2 a source size of 7.6 " (Remijan et al. 2003) is 
applied and results in Atot = 5.8 x 10'® ctn“^, and Trot = 140 K 
for HCOOCH 3 . As in the case for W51/e2, the estimated un¬ 
certainty of the column density is 30 - 40% and 40 K for the 
temperature. Results from both sources are given in Table 3 and 


we have checked these results by using the parameters derived 
from the rotation diagram. 

3.3. (CH20H)2 

Numerous (CH 20 H )2 lines are blended with other species or 
have line intensities below 3 cr. A total of two and five reason¬ 
able well separated lines above the 3 cr limit toward G34.3+0.2 
and W51/e2 respectively were chosen for this analysis (see Ta¬ 
ble 2 and 1 for the spectroscopic data and Figure 1 and 2 for 
the zoom-in of the spectra around the lines). With only 2-5 un¬ 
blended lines, the assignment of (CH 20 F [)2 cannot be consid¬ 
ered as a firm detection. However, even if our tentative detection 
of (CH 20 H )2 is not confirmed, this measurement represents a 
useful upper limit for the column density of (CH 20 H) 2 , which 
still can provide important information when compared to the 
HCOOCH 3 detection. 

When only a few data points are available, a statistically reli¬ 
able result of the rotational temperature cannot be obtained from 
the traditional use of the rotational diagram method. Thus, we as¬ 
sume that (CH 20 H )2 and HCOOCH 3 emit at the same Tex which 
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Fig. 2: Same as for Figure 1 but for G34.3+0.2. The transitions are located at 240807.88 and 242656.23 MHz - marked by black vertical lines in 
the plot. 


in LTE is equal to Tjof The uncertainty of the column density is, 
as for HCOOCH 3 , estimated to 30^0%. 

Using ^source = 7.6" and Trot = 140 K, a column density 
of 1.9 X 10'^ cm ~2 derived for (CH 20 H )2 toward G34.3+0.2. 
For W51/e2, a source size of 2.4" and Trot = 120 K returns 
Ntot = 3.1 X 10*® cm“^ for (CH 20 H) 2 . These results, listed in Ta¬ 
ble 3, are used as input parameters for synthetic spectra in order 
to perform a sanity check. We checked each line in the synthetic 
spectra against the observed spectra for each source and they 
all seem to provide a reasonable match, at least within the esti¬ 
mated uncertainty. Figure 1 and 2 shows the comparison of the 
synthetic spectra against the observed spectra. It is evident from 
Figure 1 that the synthetic spectrum in the 1 mm data toward 
W51/e2 slightly overproduces the observed spectrum, while it 
reproduces the observed line in the 3 mm data. As this line has a 
lower Fup this could indicate that a lower excitation temperature 
would give a better ht. Indeed a hxed Trot = 70 K in the rotational 
diagram returns the same column density, Ntot = 3.1 x 10*®cm“^. 

3.4. CH 2 OHCHO 

Several lines in the data might be assigned to CH 2 OHCHO. 
However, they are either blended with the emission of another 
molecules or too faint (i.e. below the 3 cr limit), thus we cannot 
claim a detection. Nevertheless, we are able to obtain an upper 
limit for the column density. Following the assumptions stated 
in section 3.1, synthetic spectra were generated using ^source, Trot 
and Av in Table 3, allowing Ntot to vary. More specifically, the 
upper limit of the column density is then determined by increas¬ 
ing the value until the synthetic spectra would overproduce the 
observed spectra intensities at the locations of the CH 2 OHCHO 
lines in the observed spectra. The upper limit on the column den¬ 
sity is < 3 X 10*"*cm“^ toward G34.3+0.2 and < 2 x 10*®cm“^ 
toward W51/e2. The uncertainty of the column density is, as for 
HCOOCH 3 , estimated to 30-40%. 

3.5. Relative abundances 

Following the three assumptions listed in in section 3.1, it 
is possible to calculate the relative abundance of the species 
at the rotational temperature and source size found for each 
source. The following abundance ratios have been com¬ 
puted; HC 00 CH 3 /(CH 20 H) 2 , HCOOCH 3 /CH 2 OHCHO and 
(CH 20 H) 2 / CH 2 OHCHO. As CH 2 OHCHO is not detected, it 
is only possible to set an upper limit on the column density. 
The HCOOCH 3 /CH 2 OHCHO and (CH 20 H) 2 / CH 2 OHCHO ra¬ 


tios are therefore lower limits. All three relative abundance ratios 
for W51/e2 and G34.3+0.2 are listed in Table 4 along with pre¬ 
vious measurements toward high-mass stars-forming regions, a 
hot core, an intermediate-mass protostar, low-mass protostars, 
molecular clouds toward the Galactic Centre and comets. 

4. Discussion 

(CH 20 H )2 has previously been detected in several other sources 
(Sgr B2(N), NGC 7129 FIRS2, NGC 1333 IRAS 2A and Orion- 
KL), with additional marginal/tentative detections in W51 el/e2 
and IRAS 16293-2422 (Hollis et al. 2002; Fuente et al. 2014; 
Maury et al. 2014; Coutens et al. 2015; Kalenskii & Johans¬ 
son 2010; Jprgensen et al. 2012). In this paper we have iden¬ 
tified (CH 20 H )2 and HCOOCH 3 in two high-mass protostel- 
lar sources, W51/2 and G34.3+0.2, and in addition been able 
to give an upper limit for the column density of CH 2 OHCHO 
in the same sources. The results from this study are roughly 
similar to the estimates ratios seen in Orion KL (Favre et al. 
2011; Brouillet et al. 2015) for the HCOOCH 3 /CH 2 OHCHO 
and HC 00 CH 3 /(CH 20 H )2 ratio (table 4) while the upper lim¬ 
its for (CH 20 H) 2 /CH 20 HCH 0 are similar to the results for 
NGC 1333 IRAS 2A (Coutens et al. 2015) as well as the upper 
limits results for comets (Biver et al. 2014; Bockelee-Morvan 
et al. 2000; Crovisier et al. 2004a,b). 

When investigating the conditions leading to differences in 
observed abundance ratios, it is important to take both the forma¬ 
tion process as well as the destruction processes of the molecules 
into consideration. Garrod et al. (2008) combine a gas-grain 
chemical network with a physical model to test if the chemistry 
can reproduce the observed abundances of previously detected 
organic molecules in physical conditions characteristic to star¬ 
forming regions. The physical model used in Garrod et al. (2008) 
is based on Viti et al. (2004) and consist of an isothermal collapse 
followed by a warm-up phase with temperatures from 10 K to 
200 K, assuming absolute timescales for the warm-up phases are 
1X10®, 2x 10® and 5 x 1 O'* yr representing low-, intermediate- and 
high-mass star-formation respectively. However Aikawa et al. 
(2008) argue that the relation should be reversed, as the warm¬ 
up timescale should be relative and should depend on the ratio 
of the size of the warm region to the infall speed instead of the 
overall speed of star formation. Either way, Garrod et al. (2008) 
and Aikawa et al. (2008) agree that the timescales, whether abso¬ 
lute or relative, at the different temperature-ranges are important 
for the chemistry. 

For all the sources in Table 4, HCOOCH 3 is more abun¬ 
dant than CH 2 OHCHO, from a factor of > 550 in W51/e2 
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Table 4: Relative abundances of (CH 20 H) 2 , HCOOCH 3 and CH 2 OHCHO in different sources 


source 

HC 00 CH 3 /(CH 20 H )2 

HCOOCH 3 /CH 2 OHCHO 

(CH 20 H) 2 / CH 2 OHCHO 

^bol[^ 0 ] 

References 

W5 l/e2 (this study) 

35 

>550 

> 16 

4.7 X 10'’ 

1,2 

Orion-KL^'') 

35 

> 200-300 

> 12 

10 ^ 

3,4 

G34.3-H0.2 (this study) 

31 

> 193 

> 6 

2.8 X 10^ 

1,5 

Sgr B2(N) LMH® 

- 

52 

- 

10 ’ 

6,7 

G31.41 -hO.31® 

- 

< 34 

- 

1.8 X 10^ 

8,5 

NGC 7129 FIRS2® 

~ 15 

~ 30 

~2 

500 

9 

IRAS NGC 1333 2A® 

4 

20 

5 

20 

10,11 

IRAS 16293-2422('') 

~ 13 

~ 13 

~ 1 

27 

12,13 

IRAS NGC 1333 4A® 

- 

10 

- 

7.7 

14 

MC G -H 0.693(^> 

4.3 

5.2 

1.2 

_ 

15 

MCG-0.11® 

2.8 

4.3 

1.6 

- 

15 

MC G - 0.02® 

2.5 

3.3 

1.3 

- 

15 

Lemmon® 

< 0.7 

- 

> 3 

- 

16 

LoveJoy^^' 

< 0.6 

- 

> 5 


16 

Hale-Bopp*^' 

0.32 

> 2 

>6.25 

- 

17 


Notes. The table show the results from similar studies of : The Orion Kleinmann-Low nebula, a high-mass star-forming region *** the Large 
Molecule Heimat (LMH) hot core source toward Sgr B2(N), ^‘^'high-mass star-forming region, intermediate-mass protostar , low-mass Class 
0 protostars, molecular clouds in the Sgr A complex (MC G - 0.02 - 0.07 ("the 50 km s“' cloud") and MC G - 0.11 - 0.08 ("the 20 km s“' 
cloud")) and in the Sgr B2 complex (MC G -t- 0.693 - 0.03) and comets. 

References. (1) this study; (2)Urquhart et al. (2014); (3) Favre et al. (2011); Brouillet et al. (2015); (4) Crockett et al. (2014); (5) van Dishoeck 
et al. (2011); ( 6 ) Hollis et al. (2001); (7) Goldsmith et al. (1990); ( 8 ) Beltran et al. (2009); (9) Fuente et al. (2014);(10) Coutens et al. (2015); 
(11) Jprgensen et al. (2009); (12) Jprgensen et al. (2012, in prep.); (13) Schoier et al. (2002); (14) Taquet et al. (2015); (15) Requena-Torres et al. 
(2008); (16) Biver et al. (2014); (17) Bockelee-Morvan et al. (2000); Crovisier et al. (2004a,b). 


to a factor of > 2 in comet Hale-Bopp. According to Garrod 
et al. (2008) HCOOCH 3 and CH 2 OHCHO have similar forma¬ 
tion pathways, which are based on the addition of HCO and 
CH 3 O or CH 2 OH at 30^0 K. Garrod et al. (2008) find that 
the production rates of CH 3 O and CH 2 OH are the same, which 
leads to similar abundances for HCOOCH 3 and CH 2 OHCHO. 
Intuitively this makes sense as the two molecules are isomers 
but it contradicts observations as HCOOCH 3 is observed to be 
much more abundant than CH 2 OHCHO in all the sources re¬ 
ported so far. As suggested by Garrod et al. (2008), this dis¬ 
crepancy could be due to the differences in the CH 3 O/CH 2 OH 
branching ratio. According to Oberg et al. (2009), HCOOCH 3 
forms at a lower temperature than CH 2 OHCHO, which can also 
have a impact on the resulting ratio of the two species. An¬ 
other explanation to the large HCOOCH 3 /CH 2 OHCHO ratio is 
the assumptions regarding thermal evaporation for these two 
species. Garrod et al. (2008) suggest that CH 2 OHCHO remains 
on the grains until is co-desorbs with water (at ~ 110 K) while 
HCOOCH 3 evaporates at 70-80 K. This leaves CH 2 OHCHO to 
be destroyed by OH radicals at higher temperatures on the grains 
prior to evaporation. As Oberg et al. (2009) conclude, the ratio of 
HCOOCH 3 /CH 2 OHCHO does not depend greatly on the initial 
ice composition, and so the observed variations will most likely 
be linked to the different temperature conditions of the different 
sources. 

In Table 4 the sources are listed in order with decreas¬ 
ing HCOOCH 3 /CH 2 OHCHO ratio, which in turn also roughly 
correspond to a descending order of luminosity, which is 
also listed in the table. The top plot in Figure 3 shows a 
schematic bar-plot of HCOOCH 3 /CH 2 OHCHO against Lboi- 
On the x-axis the sources in Table 4 have been plotted in 
descending order from left to right, but not to scale. As 
illustrated by Figure 3 (top) a rough correlation between 
HCOOCH 3 /CH 2 OHCHO and Lboi exist: more luminous sources 


show a high HCOOCH 3 /CH 2 OHCHO ratio, while low lumi¬ 
nosity sources show a low HCOOCH 3 /CH 2 OHCHO ratio with 
intermediate values in between. The apparent correlation be¬ 
tween HCOOCH 3 /CH 2 OHCHO and source luminosity supports 
the hypothesis that the HCOOCH 3 /CH 2 OHCHO ratio depends 
mainly on the temperature, which in turns depends on the lumi¬ 
nosity of the source. Even if the chemistry and the temperature 
profile in all the sources were quite similar, the differences in 
timescales at the different temperature-ranges (low, intermediate 
to high temperature) are not. A higher HCOOCH 3 /CH 2 OHCHO 
ratio in more luminous sources, could simply be a consequence 
of i) more luminous sources have experienced a longer timescale 
at temperatures which are either more favorable to the forma¬ 
tion of HCOOCH 3 and/or to the destruction of CH 2 OHCHO 
or ii) less luminous sources experience a shorter timescale at 
high temperature, which would conserve a greater fraction of 
CH 2 OHCHO than their more luminous counterparts. One should 
of course keep in mind that the sources discussed here have tem¬ 
perature profiles which covers a range of temperatures falling of 
with radius and that hot core regions are permeated by shocks 
which also produce a range of temperatures. For the Galac¬ 
tic Centre molecular clouds from the study by Requena-Torres 
et al. (2008) no luminosities are listed: these are warm, low den¬ 
sity clouds with no sign of star-formation. However, judging by 
the HCOOCH 3 /CH 2 OHCHO and HC 00 CH 3 /(CH 20 H )2 ratios, 
they appear to be closer to the low-mass protostars than to the 
hot cores and high-mass star-forming regions. 

In the bottom plot in Figure 3, the HC 00 CH 3 /(CH 20 H )2 
abundance ratio also shows a decrease with Lboi similar to that of 
HCOOCH 3 /CH 2 OHCHO. If the same temperature timescale ar¬ 
gument is applied to this correlation, one would expect an oppo¬ 
site trend, as (CH 20 H )2 is formed a high temperatures (Garrod 
et al. 2008). A possible explanation for this seemingly contradic¬ 
tion is that less luminous sources might experience high temper- 
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Fig. 3: Schematic bar-plot of HCOOCH 3 /CH 2 OHCHO (top) and 
HC 00 CH 3 /(CH 20 H )2 (bottom) against L(,oi- The sources from Table 4 
have been plotted in descending order of luminosity from left to right 
on the x-axis, but not to scale, and the sources have been grouped into 
high-mass sources (blue), intermediate-mass protostar (red) and 
low-mass protostars (green). The white bars for four of the high-mass 
sources are upper/lower levels, which is indicated by the direction of 
the arrow. The luminosity for the top plot spans from Sgr B2(N) with 
Lboi = lO^L© to IRAS NGC 1333 4A with Lboi = 7.7 Lq, while it spans 
from W51/e2 with Lboi = 4.7 x 10®L 0 to IRAS NGC 1333 4A with 
Lboi = 20Lq in the bottom plot. 


atures at timescales, that are just long enough for (CH 20 H )2 to 
form, but not long enough for (CH 20 H )2 to be destroyed again. 
Another possible explanation is that the initial ice composition 
affects the abundance ratio. While the HCOOCH 3 /CH 2 OHCHO 
ratio is independent of initial ice composition, then the formation 
of (CH 20 H )2 is strongly correlated to the CH 30 H:CO composi¬ 
tion of the ice (Oberg et al. 2009). Oberg et al. (2009) show that 
pure CH 3 OH ice strongly enhances the (CH 20 H )2 abundance 
as compared to ice mixes containing CO. This contradicts the 
predictions by Garrod et al. (2008) where the (CH 20 H )2 abun¬ 
dance actually drops 1-2 magnitudes when the initial CH 3 OH 
ice-composition in their model is reduced by a factor of 10. How¬ 
ever, to date there is no observational evidence that the ice con¬ 
tents on the grains vary statistically significant from high-mass 
protostars to low-mass protostars (Oberg et al. 2011). 

Until recently a high (CH 20 H) 2 /CH 20 HCH 0 ratio was 
speculated to be specific property of comets (Biver et al. 2014). 
But as Coutens et al. (2015) has reported a values of ~ 5 for a 
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low-mass protostar and we, together with Brouillet et al. (2015), 
report lower limits of > 6-16 for high-mass sources, then a 
(CH 20 H) 2 /CH 20 HCH 0 ratio larger than 3 can be expected to 
be observed in other sources as well. 


5. Conclusions 

In summary, we have tentatively detected (CH 20 H )2 in 
G34.3H-0.2 for the hrst time and our tentative detection in 
W51/e2 confirms the previous marginally detection by Kalen- 
skii & Johansson (2010). In addition, we derive upper lim¬ 
its for the column density of CH 2 OHCHO emission in both 
sources. From these data we determine the relative abundances 
of (CH 20 H) 2 , X and CH 2 OHCHO. The relative abundances 
of these species are compared to measurements from litera¬ 
ture covering a wide range of source environments and lumi¬ 
nosities. The data show what appears to be a correlation be¬ 
tween source luminosity and HC 00 CH 3 /(CH 20 H )2 as well 
as HCOOCH 3 /CH 2 OHCHO. This apparent correlation is pro¬ 
posed to be a consequence of the relative timescales each source 
spend at different temperatures-ranges in their evolution. Using 
the upper limit for the column density of CH 2 OHCHO gives 
a lower limit for (CH 20 H) 2 /CH 20 HCH 0 of > 16 and > 6 
for W51/2 and G34.3H-0.2 respectively. These results, together 
with an upper limit of > 12 for Orion-KL (Brouillet et al. 
2015) and (CH 20 H) 2 /CH 20 HCH 0 = 5 for IRAS NGC 1333 
2A (Coutens et al. 2015), shows that one can expect to find high 
(CH 20 H) 2 /CH 20 HCH 0 abundance ratios in multiple types of 
environments. 

Additional systematic surveys of these and other relevant 
molecules in additional sources as well as more model and lab¬ 
oratory work is needed in order to fully constrain the formation 
pathway of complex molecules. In particular, the Atacama Large 
Millimeter/submillimeter Array (ALMA) shows great potential 
for successfully revealing the formation processes with its high 
sensitivity and resolution making it possible to map the relative 
spatial distributions of these sources. 
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Appendix A: Reference model 


Table A.l: Column densities of species in our reference model 


Molecule 

G34.3+0.2 

W51/e2 


[ctn“^] 

[ctn“^] 

CH 3 OCH 3 

4x 10 '® 

3 X 10'^ 

HCOOH t 

1 X 10'® 

1 X 10'® 

CH 3 CHO 

2 x 10 '® 

1 X 10'® 

CH 3 OH 

1 X 10'^ 

3x 10 '*' 

C 2 H 5 OH 

1 X 10'® 

1 X 10'^ 

CH 3 CN 

2 x 10 '® 

6 x 10 '® 


Notes. Estimated column densities of other relevant organic species. 
These values are used to create synthetic spectra of these molecules, 
which serves as a reference model when identifying lines. 
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Fig. B.l: Synthetic spectra (red: (CH20H)2, yellow: HCOOCFI 3 , green: CH 2 OHCHO, blue: the molecules listed in Table A.l) on 
top of observed spectrum (black) at 1 mm for W51/e2 in the frequency range ~238-242 GHz. 
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Fig. B.2: Same as Figure B.l but for ~242-246 GFIz. 
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Fig. B.3: Synthetic spectra (red: (CH 20 H) 2 , yellow: HCOOCH 3 , green: CH 2 OHCHO, blue: the molecules listed in Table A.l) on 
top of observed spectrum (black) at 1 mm for G34.3+0.2 in the frequency range ~238-242 GHz. 
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Fig. B.4: Same as for Figure B.3 but for ~242-246 GHz. 
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Fig. B.5: Synthetic spectra (red: (CH20H)2, yellow: HCOOCH 3 , green: CH 2 OHCHO, blue: the molecules listed in Table A.l) on 
top of observed spectrum (black) at 3 mm for W51/e2 in the frequency range ~100-103 GHz. 
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Fig. B.6: Same as for Figure B.5 but for ~103-106 GHz. 
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Fig. B.7: Synthetic spectra (red: (CH 20 H) 2 , yellow: HCOOCH 3 , green: CH 2 OHCHO, blue: the molecules listed in Table A.l) on 
top of observed spectrum (black) at 3 mm for G34.3+0.2 in the frequency range ~100-103 GHz. 
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Fig. B.8: Same as for Figure B.7 but for ~103-106 GHz. 
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